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an external field coupled to the last matrix of the chain) . We show that the topological expansion of the free 
energy, is, like for the 1 and 2-matrix model, given by the symplectic invariants of |19| . As a consequence, we 
find the double scaling limit explicitly, and we discuss modular properties, large N asymptotics. We also briefly 
discuss the limit of an infinite chain of matrices (matrix quantum mechanics) . 



tlnstitut de Physique Theorique, CEA, IPhT, F-91191 Gif-sur-Yvette, France, 
CNRS, URA 2306, F-91191 Gif-sur-Yvette, France. 

X Centre de Recherches Mathematiques, Pavilion Andri'JEj-Aisenstadt, Universiti'JEj de Montri'JCjal. 
2920, Chemin de la tour, Montrf£jal (QuidEjbec) H3T 1J4, CANADA. 

1 E-mail: bertrand.eynard@cea.fr 

2 E-mail: pratsferrer@crm.umontreal.ca 



1 



1 Introduction 



Since the famous discovery of Brezin, Itzykson, Parisi and Zuber [7], it has been known and widely used, that 
formal matrix integrals are generating functions for the enumeration of discrete surfaces of given topologies (the 
role of topology was first noticed by 't Hooft [35] )■ The 1-matrix model is known to count discrete surfaces 
obtained by gluing polygonal pieces side by side. It is the partition function of random discrete surfaces [13], 
also called random "maps". 

Other matrix models are also partition functions of random discrete surfaces, with additional "colors" on 
the faces [251 [13]. 

In particular, the "2-matrix model" is a partition function of random discrete surfaces, whose polygonal 
pieces can have two possible colors (or say, two possible spins + or -), and surfaces are counted according 
to the number of edges separating polygons of different colors, that is polygons with different spins. Thus it 
counts surfaces with a weight proportional to the exponential of J2<ij> a i (J j (where the sum is over pairs of 
neighboring pieces, and Oi is the spin of the piece i). In other words this is an Ising model on a random discrete 
surface [28] . 

The most natural generalization is the "Chain of matrices" matrix model. It is the generating function for 
counting discrete surfaces, where pieces can have a color i € [1, . . . , n], and where each spin configuration on 
the surface is weighted by Yi<i j>(C~ 1 )i,j where C is a Toeplitz matrix of the form: 
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( & 

-Cl,2 



-Cl,2 

(2) 
92 



"C2,3 



o \ 



C-n — 1 , n 

9i n) J 



(1.1) 



\ -Cn-l, r , 

The partition function for the chain of matrices is the formal small T expansion of the following matrix integral: 

(1.2) 
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where M g (v) is the set of all connected orientable discrete surfaces of genus g with v vertices, with m,k polygonal 
pieces of size k (i.e. k— angles) of color i, and n e dges<i,j> edges separating colors i and j, and #Aut is the number 
of automorphisms of the surface. Notice that for fixed g and v, M. g (v) is a finite set, and therefore F g is indeed 
a formal series in T. 

One may also be interested in discrete surfaces with m marked faces, whose generating function is given by: 
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where M^^...^ is the set of all connected discrete surfaces of genus g, with n^k A:— angles of color i, and 
fiedges<ij> edges separating colors i and j, and with m marked faces (and with one marked edge on each 
marked face), of respective perimeters h, . . ■ ,l m , and colors ii, ■ . ■ ,i m . Again, for fixed m and g, there are 
finitely many such surfaces with a given number of vertices, and the sum is a formal power series in T. Notice 
that if there is only one marked face m = 1, i.e. one marked edge, we have a rooted map, and #Aut(5) = 1. 

Recently, the computation of the F g 's and W^'s was completed for the 1-matrix model (n = 1) in [161 18]. 
and 2-matrix model (n = 2) [9J [THl [19] , and our goal is to extend the method of [H] to the chain of matrices of 
arbitrary (but finite) length n > 1. 

In fact, the method of [19] allows to find the solution for a generalization of the chain of matrices, where in 
addition, the last matrix is coupled to a fixed matrix M n+ i, called external field. Matrix models with external 
fields also have some combinatorial interpretations, and have been studied for various applications. The most 
famous is the Kontsevich integral, which is the generating function for intersection numbers [33l [24] fl9] . 

Here, we solve this more general model. 

Multimatrix model also play an important role in quantum gravity and string theory, where they play the 
role of a regularized discrete space-time. The 1-matrix model, counts discrete surfaces without color, and is a 
model for quantum gravity without matter, whereas the chain of matrices counts discrete surfaces with n colors, 
and is interpreted as a model of quantum gravity with some matter field |27 [ fT0 ] fT][29 ] [30 l ll3]. namely a matter 
which can have n possible states. More recently, matrix models have played a role in topological string theory 



Outline of the article: 

• In section[2]we introduce all the definitions and notations necessary for the derivation of the loop equations. 
These are quite clearly inspired by the work on [15] where the loop equations were already found in a 
slightly less compact way. 

• In section [3] we derive the master loop equation that will allow us to solve the model. We also consider 
the -^2 expansion here and find the spectral curve for this model. 

• In section |4] we overview all the important algebraic geometry tools and the algebraic curve properties 
that are relevant for us. 

• In section [5]we apply the same techniques of [9] to prove uniqueness of the solution and to find the actual 
solution for the correlators of the first matrix Mi of the chain. 

• In section [6] we find the variation of the curve, and all the correlation functions, in terms of the moduli 
of the chain of matrices. This leads us to an expression for the whole topological expansion of the free 
energy for the chain of matrices. 

• In section \7\ we study some corollaries of the properties of the symplectic invariants of [19] , in particular 
we get the double scaling limit, and modular properties. 

• In section [U we briefly discuss the "matrix quantum mechanics", i.e. the limit of an infinite chain of 
matrices. 

• Finally, section [9] is the conclusion. 
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2 Notations and Definitions 



2.1 The formal chain matrix model with external field 

The formal chain matrix model with external field, is a formal matrix integral!, with n matrices of size N with 
potentials Vj(M,-), arranged in a chain with Itzykson-Zuber like interactions: 

/n 
J] dM z e-H^Wl-wW) (2.1) 
«— i 

where M n +i is a constant matrix, which we may choose diagonal M n +i = A , with s different eigenvalues A, 
and multiplicities k h = N): 

A = diag( Ai, ■ ■ ■ , Ai , • ■ • , A;, • ■ ■ , A,, ■ • ■ , A g , •■ • , A s ). (2.2) 
h U h 

It reduces to the usual "Chain of Matrices" when A = M n+ \ = 0. 

The measures dMj = YIa=i ^-^jj Y[^<k ^(M^ )d^s(Mjl! ) are the usual Lebesgue measures for hermitian 
matrices. The potentials Vi(x) are polynomials^ of degree di + 1, 

di+l (i) 

v ^ x ) = E 9 ~T xk ( 2 - 3 ) 

fc=l 

but the same results contained in this paper can clearly be extended to functions Vi whose derivatives V{ are 
rational functions. In general we are interested in formal expectation values of functions of Mj defined by 

(/(Mi, • • • , M n )) = _L /" TT dMi /(Mi, • • • , M n ) e-^H^Uiv^-c^M.M^)) (2 _ 4) 
Z Ch J f \ 

but we will also be interested in the free energy defined as the logarithm of the partition function Zch- 

The expansion can be considered when we work with the formal version of this matrix integral. What 
that means is that we must interpret the integrals as a formal expansion of all the non-gaussian terms in the 
exponential and perform the integral as a perturbation integral around a minimum of the action 

tr Vi{Mi) - Ci.i+iMM+iJ . (2.5) 

The equations that define a minimum are 

Vi'(Mi) = ci, 2 M 2 , T4'(M fe ) = Cfc_i, fc Mfc_i + c k ,k+iM k +i k>2 (2.6) 
In particular we can choose a minimum such that all Mfe's are diagonal Mk = diag(/2i , . . . ,p, N ), which satisfy: 



nWi j — cififJ-i 

jr'fr( k )\ -C=-l) , -0+1) , r, 

VkKN ) = Ck-i,kN + Cfe,fc+iMi fe = 2, • • • ,n 



(2.7) 



with = = Note that c„.„+i can be absorbed into A, so that we will fix it to 1. These equations 

have D = d\d,2 ■ ■ ■ d n s solutions. Choosing which minimum we are going to perturb around, means choosing 
how many eigenvalues we are going to put on each of the D different solutions. Let us call these rii, • • • , no, 
with the restriction ^ ■ n, = N. In the following we are going to refer to e, = T ^ as the filling fractions. 

In other words, for each choice of filling fractions e = (ei, . . . , cd-i), we can define a formal integral by 
perturbation around the corresponding minimum. Almost by definition, there must exist anti-clockwise contours 
Ai, i = 1, . . . , D, such that 

7 tr ( —!—) \ dx = a =T^ (2.8) 



2inN J A . \ \x-MxJ J N 



3 A formal integral is defined as the exchange of the integral and the Taylor expansion of the exponential of non-quadratic terms 
in the potentials, see |23j . 

4 Notice that here, in contrast to equation 111. 3(1 . we allow for a linear term in the potential. This is convenient but can be 
trivially undone by a shifts proportional to the identity. 
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2.2 Definitions of correlation functions 



In order to define the good observables of our model, we first need to introduce (like in [15]) the following 
polynomials fi,j(x\, . . . , x n ) 



fi,j( x i, > X j) — J 



1 



k — i 



Ck-l,k 



■det 



V!(xi) -Cis + ix i+ i 
-Ci,i+iXi V{ +1 {x i+ i) 



-Cj.j+iXj-i V'(xj) 



if i < j 



(2.9) 



1 if' j I I 

= if i > j + 1 

They satisfy the recursion relation 

Ci-l,ifi,j(Xi, . . . ,Xj) = Vl{Xi)fi + ij(Xi + u . . . ,Xj) - Ci t i +1 XiXi + ifi + 2,j(Xi + 2, ■ ■ ■ , Xj) 



(2.10) 



with the initial conditions fk+i,k = 1, and fk+i,k = for all I > 1. The first polynomials generated by this 
recursion relation are 



fi,i [Xi) 



*7(si) 



/i— 1,2 V^i— 1 5 ) — Xi — lXi 



fi-2,i( x i-2: Xi-i, Xi) 



Define also the following functions 



<H-2,i-l c i-l,i Ci-2,i-l 

VU^i-2) V?-i(ai-i)V7(si) 

Cj-3,i-2 Cj_2,i-1 Cj-i^ 
^-2(^-2) Ci_i,i _ Cj_2,»-1 Kfai) 



(2.11) 



Cj-3.i-2 



Ci-1, 



U>i(xi) = 



1 



- Ah 

where S(z) is the minimal polynomial of A: 



Q(z) 



1 S(z) - S(A) 

Cn,„+1 Z — A 



S(z) = l[(z-X i ) 



(2.12) 



(2.13) 



4 



The loop equations in following sections will be written in terms of the following matrix model observables or 
correlation functions. 

Wo(xi) = ^tr(wi(a:i)) 
P(xi) = Pol/i,i(xi)W (a:i) = Po\V(( Xl )W (xi) 

X\ X\ 

W i (x 1 ,x i , ...,x n ,z)= Pol fiJxi, . . . ,x n ) ( -77 tr (wi(xi)wi(xi) ■ ■ ■ w n (x n )Q(z)) ) , 
for i = 2, . . . ,n — 1 

/ T 

Wx(xi,...,a:„,z) = Pol /i, n (xi, . . . ,£„) ( — tr (u>i(a;i) ■ • ■UMZnjQO)) 

Xl,...,x n \ iV 



(2.14) 



— > is a polynomial in all variables 
/ T 

W n (xi,z) = ( — tr (wi(xi)Q(z)) 

W 0; i(xi;a?i) - dy^ WoOEi) = (tr (wi(xi)) tr (w 1 (x' 1 ))) c 
d 

Wi-llXi , Xi , . . . , X n , Z\ Xi) = . rrWi (x\ , Xi , . . . , X n , z) = 

oVi{x[) 

tr^i^i)) Pol fi n (xi,...,x n )tr(wx(x 1 )wi(x i )---w n (x rl )Q(z)) 

Xi,...,X n 

where the symbol ll Pol x f(x)" represents the polynomial part on the variable x in the vicinity of 00 of the 
function "f(x)", and the loop insertion operator simbol is defined by 

d Id y-^ k d 

1) ~ 1^ Zfc+TTTIT ( 2 - 15 ) 



9Vi(x) x dg (f) ^x^ d £) 
like in [15]. At some point we will write the topological expansion^! of some of these functions, for example 



00 / rp \ 2h 

W ^) = J2(n) W o H) M (2.16) 

h=0 ^ ' 

and similarly for other functions. These are all the definitions we need for the derivation of the loop equations 

3 Master Loop Equation 

To find the master loop equation (proceeding as in [15]) we are going to consider the following local changes of 
variables 

SMi = e Pol fi+i, n (x i+1 , . . . ,x n )w i+ i(x i+1 ) ■ ■ ■ w n {x n )Q(z)w\{xx) + 0(e 2 ), l<i<n 

Xi+l,...,Xn ^ 

6AI n = eQ{z)w 1 (x 1 ) + 0{e 2 ) 

with e a small parameter. Notice that 5Mi does not contain Mi except for i = 1. We must then consider 8M\ 
separately. 

3.1 Loop Equation for 8 Mi 

Consider the change of variables 

5Ah=e Pol f2,n( x 2,---,x n )w2{x 2 )---w n (x n )Q(z)w 1 (x 1 ) + 0(e 2 ) (3.2) 

X2,...,X n 



6 The topological expansion of a formal integral, is not a large TV expansion, it is a small T expansion, and for each power of T, 
the coefficient is a polynomial in iV~ 2 . The Wq H \x) is merely the formal series in T, containing the degree h terms. 
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The first order variation in e of the integral 1|2.1| gives the Schwinger-Dyson equation (called loop equation in 
the matrix model context): 

(j^2 tT (Wi(xi))tr \ Wl{xi) ^ Pol^ f2,n{x2, ■ ■ ■ ,X n )w 2 (x 2 ) ■ ■ ■ W r Xx„)Q{z)^j 

= \T7 tv (wi(xi) (V({Mi) ~ ci, 2 M 2 ) Pol fa, n (x2, . . . ,x n )w 2 (x 2 ) ■ ■ ■ w n {x n )Q(z) 

\JS \ x 2 ,...,x rL 



(3.4) 



(3.3) 

JS n II 

Using (|2.10p we find, after some algebra, the loop equation 
T 2 

j^W 2 -i{x\,x 2 , ...,x n ,z; x\) + [01,2x2 - V{(xi) + Wo(x\))W 2 (xi, x 2 , ...,x n ,z) = 

= - Wi(xi, . . . ,X n , z) + (V 2 (x 2 ) ~ Ci, 2 Xi)W 3 (xi,X 3l . . . ,X n , z) - C 2 ^X 2 Wi(x\,XA, ...,X n ,z) 
- (j^ tT f^Xl) (Va(M 2 ) - Cl, 2 Ml) ^Pol_ f 3 , n {x 3 , ■ ■ ■ ,X n )w 3 (x 3 ) ■ ■ ■ W n (x n )Q(z)^ 

+ \TFtr (wi [x 1 )c 2)3 M 2 Pol / 4 ,„(x 4 , . . . , x„)w4(x 4 ) ■ ■ • w n {x n )Q{z) 

\JS \ Xi,...,X n 

3.2 Loop Equation for 5 Mi 

The rest of the loop equations follow the same principle. We will compute the remaining in one shot. 

SMi = e Pol f i+lt , l (x l+1 ,...,x n )w l+1 (x l+1 )---w n (x n )Q(z)w 1 {x 1 ) + 0(e 2 ) (3.5) 

from which the order e variation of the partition function is 
= (^tr (wi(xi) (V-(Mi) - Cj_i,iMj_i) ^ Pol_ f i+1} „(x i+ i, . . . , x„)w i+1 (x i+ i) ■ ■ ■ w„(x n )Q(z)J ^> 

- Ci >i+ xXi+xWi+l{xx,Xi+l, ... ,X„,z) + V- +1 (x i+ i)Wi + 2{x 1 ,Xi + 2, ...,X n ,z) — Ci,i+lXi+lWi+3(xl,Xi+3, ■ ■ ■ ,X„,z) 

I T ( \\ (3.6) 

- ( ]y tr [ Wl ( Xl ) V i ( Mi ) Po1 fi+2,n{x l+ 2, ■ ■ ■ , X n ) W l+2 (X i+ 2 ) • • • W„(x n )Q(z) j ) 

+ (jjtr \ U)l(Xl)Ci,i + lM i+ l ^ Pol ^ fi+3,n(xi+3, ■ ■ ■ , X n )w i+3 (x i+3 ) ■ • ■ W n (x n )Q(z)^j \ 

In particular, for i = n due to the fact that f n +2,n = fn+ 3 .n — we have 

= ^tr (wi(xi) (K(Mn) - Cu-lA-i) Q(z))\ 
zW n+1 (x 1 ,z) + S(z)W (x 1 ) 

3.3 Master Loop Equation 

When we sum up equations l|3.4p and l|3.6p for i = 2, . . . , n we find the master loop equation 

-r75-W 2 ;l(xi,X2, . . . , X„, X n+ i; Xl) + (a,2X 2 - V{(xi) + W (xi)) (W 2 (xi,X2, ■ ■ ■ ,X„,X n+1 ) - S(x n +l)) 

= - Wl(xi, . . . ,X n ,X n +l) + (V[(xi) - 01,2X2) S(x n +l) 
n 

+ y)(yi(xj) — <k-i t iXi-i — Ci ii+ iXi + i)Wi + i(xi,xi + i, . . . ,x„+i) 



(3.7) 



(3.8) 



where we have redefined z = x n+ \. Remember that W\{x\, • • • , x n ) is a polynomial in all its variables and that 
) is a polynomial in all its variables except x\. In particular, we may choose Xi = x i (x ll x 2 ), 
i = 3, ... ,11 such that 

V-(xi) = Ci-ijXi-i + c hi+1 x i+1 , Vi = 2, . . . , n (3.9) 
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and in that case (|3.8p reduces to 
T 2 

— W 2 ;i(xi,x 2 ;xi) + (ci. 2 x 2 - Y(xi))U(xi,x 2 ) = - W\{x\,x 2 ) + {V((xi) - c^ 2 x 2 )S{xi,x 2 ) 

iv z (3.10) 

=E(xi,x 2 ) 

where we have defined 

Y{x 1 )=V{(x 1 )-W {x 1 ) 
U(xi,x 2 ) = W 2 (x 1 ,x 2 ,x 3 , . . .,x n +i) - S(x n+1 ) 
W 2 -i(xi,x 2 ;xi) = W 2 -i(x 1 ,x 2 ,x 3 , . . .,x n+ i;xi) (3.11) 
Wi(xi,x 2 ) = W 1 (x 1 ,x 2 ,x 3 , . . . ,x„+i) 
S(xi,x 2 ) = S(x n+1 ) 

and Xi are defined recursively from the constraints f|3.9|) 

C 2 , 3 X 3 = C2 ! 3Xs(x 1 ,X 2 ) = V 2 (x 2 ) - C la Xi 

c 3A X4 = c 3A x 4 (xi,x 2 ) = V 3 (x 3 (xi,x 2 )) -c 2 . 3 x 2 (3.12) 
Ci-i,iXi = Ci-i,iXi(x\,X2) = Vi_ 1 (x i -i(x 1 ,X2)) - Ci- 2 ,i-iXi- 2 (xx,x 2 ) , for i > 4 

Note the resemblance between equations p.9j) and l|2.7j) . 



3.4 Planar limit 

To leading order at large N, we drop the T 2 /N 2 term in the loop equation (|3. lOf) and we get: 

(c h2 x 2 - Y^( Xl ))U^{ Xl ,x 2 ) = EM( Xl ,x 2 ) (3.13) 



Notice that E^°\x\,x 2 ) is a polynomial in its 2 variables x\ and £2- 
The algebraic equation 



E^(x u x 2 ) =0 



(3.14) 



is called the spectral curve. In some sense it is the large N limit of the loop equation when we choose c\ 2 x 2 = 

y(°)(xi). 

The equation (|3.10p is of the same form as the one solved in [9] for the 2-matrix model, or the one solved 
in [19] for the 1-matrix model with external field. It can thus be solved using the same methods. Note that (as 
we said above) we consider fixed filling fractions in the formal model, which means that the 

€l = ^o 0) Mp))cMp) (3.15) 

are fixed data of the model. As a consequence of that, all the differentials WQ h \x\{p))Ax\{p) with h > 1 may 
have poles only at the branch points a>i , and all their A cycles integrals are zero due to the fixed filling fractions 
condition 

° = f A W o h) ( x ^P)) dx ^P)- ( 3 - 16 ) 

4 Algebraic geometry of the spectral curve 

The solution of the model relies on the understanding of the underlying large N spectral curve, and its algebraic- 
geometry properties. Let us see first what are the main features of E^ (xi, x 2 ) and then we will present a set 
of tools and concepts that we will need later. 

First, E^(x u x 2 ) = -W^\ Xl 

1) + (V{( x i) ~ Ci >2 x 2 )S(x n+ i) as we have noted before, is a 
polynomial in all its variables. is a polynomial of degree — 1 in the variable Xi and s — 1 in z while 

(V{(xi) — Cx <2 x 2 )S(x n +i) is clearly a polynomial of degree d\ in x±, 1 in x 2 and s in z — x n +i- 
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The relations (|3 . 1 2[) express Xi (i > 3) as a polynomial of x\ and x 2 . For example, x^{x\, X2) is a polynomial 
of degree 1 in X\ and d 2 in x 2 . In general, for i > 3, x i (x 1 ,x 2 ) a polynomial of degree YTj^3 4? in x\ and rij=2 4? 
in X2. 

With this information we see that for n > 1, (xi, x 2 ) is a polynomial of degree d\ + d 3 . . .d n s — 2 in 
x\ and a polynomial of degree d 2 ■ ■ ■ d n s — 1 in x 2 , while (V{(xi) — ci y2 x 2 )S{x n +i) is a polynomial of degree 
d\ + c?3 . . . d n s in x\ and 1 + d 2 d 3 . . . d n S in X2, i.e. 

deg Xl E m =d 1 +d 3 ...d n s = d 1 +D 1 , deg X2 £ (0) = 1 + d 2 d 3 . . . d„ s = 1 + D 2 (4.1) 

One can check from algebraic geometry usual methods (Newton's polytope for instance), that an algebraic curve 
with those degrees, has a genus ~g: 

g < d x d 2 ...d n 5 (4.2) 

So far, most of the coefficients of are not known, because they come from the unknown polynomial W\. 
However, the number of unknown coefficients of Wi, is did 2 ■ . ■ d n s — 1, and it matches precisely the generic 
genus of the spectral curve g (all the terms of Wi lie in the interior of Newton's polytope), and therefore the 
polynomial W\ (and thus E^) is entirely determined by the filling fraction conditions (we have J ^2 i q = T): 

Vi = l,...,d 1 d 2 ...d n s, €t= Jr- I Y(x)dx , E^\x, Y{x)) = (4.3) 

2j7r J A, 

Those d\d 2 . . . d n s equations determine Wi (0) and thus E^K 



4.1 Analytical structure, sheets and poles 

The algebraic curve E^'{xx,x 2 ) = has the following structure. For each value of x\ there are D 2 + 1 different 
values of x 2 , and for every value of x 2 we find D\+d\ values of x\. This observation is what defines, respectively, 
the X\ sheet structure and the x 2 sheet structure. 

The algebraic curve E^°'(xi, x 2 ) = can be parametrized as follows: there exists a compact Riemann surface 
C and two meromorphic functions x\ and x 2 on C, such that 

E(°\xi,x 2 ) = o 3p e C I xi = x 1 (p) &ndx 2 = x 2 (p) (4.4) 

Notice that the functions Xi(p) = Xi(xi(p),x 2 (p)) are also meromorphic functions on £, which satisfy: 

VpeC, V({xi{p)) = C4-\ t iXi-i{p) + Ci,i+xXi + i{p) (4.5) 

There are 5 + 1 different points on the curve where X\ (and all the other Xi(p), i < n) becomes infinite. Around 
one of these points, that we will call p = 00, a good local coordinate is (p) = x ^ p ^ ■ The 00 point in the curve is 
quite important as it marks the so called physical sheet. The other "infinity" points correspond to the situation 
where x n +\(p) = Xi, and will be called p = Xi. A good local coordinate around these points is z% (p) = x Tp) 
(a different good local coordinate could be z~ x {p) = x n +\{p) — A, which behaves as ^p—,^ - 1 Xr ^( p ) )- 
Explicitly, the negative divisor of Xk (p) is 

[x k {p)\- = -r k co - s k 2J h 

where r x = 1, r k = d x d 2 . . . d k -i 

s n +i = 0, s n = 1, Sk — d k+1 d k+2 ...d n 



Locally, near 00 we have: 



MP) = ^^~^ r]+ 0(x 1 ( P n (4.7) 
x k ( P ) ~ x^pfx (4.8) 
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and near Ai: 

p^A, C nyn+1 Jy X n+ i{p) - Ai 

1 Tl- 1 _ ' ' 

a?n+l(p) = Ai H pr + O (x„0) _ ) 

P^A, C„,„+i TV X n (p) 

a*(p) ~ (a; T[+1 (p)-Ai)- Sfc (4.10) 

P^A, 

4.2 Branchpoints and conjugated points 

From Riemann-Hurwitz, there ares+2g+si points a, on C, such that d X2 E^°'(xi, x%) = and d xl E^°'(xi, x?) ^ 
0. They are called the x\ branch points. They are the zeros of the differential dx±(p). 

For the moment, we assume that the branch points are simple, i.e. that at those points dxi(p) vanishes 
linearly when p — > a. The spectral curve is said to be regular. A spectral curve with non simple branch-points 
is called singular or critical. We study critical points below in section 17751 

Assuming that the spectral curve is regular means that near any branch-point a, Y = C\^X2 behaves locally 
like a square root Y(xi) ~ Y(x\(a)) + C yfx\ — x\{a), and therefore, for any p in the vicinity of a, there exists 
a unique point p ^ p in the same vicinity of a, such that 

x 1 (p) = x 1 (p). (4.11) 

We say that p is the conjugate point of p. The conjugate point, is defined locally near every branch-point, and 
in general it is not defined globally (see [19]) . 

4.3 Non-trivial cycles 

If C is of genus g, there exists a symplectic basis of non-trivial cycles Ai,Bj, i, j = 1, . . . ,g, such that: 

At n Bj = s hj , Ai r\Aj = o , B € n Bj = o (4.12) 

Such a basis is not unique, and we have to choose one of them. Different choices give different solutions of the 
loop equations. The choice is related to the choice of filling fractions. 

Changes of symplectic basis are called modular transformations, and, following |19| [22] we study modular 
transformations of the F g 's and Wq 8 in section [7731 

Once we have chosen a basis of non-trivial cycles, the domain £\(U;*4i U; Bt) is simply connected and is 
called the fundamental domain. 

4.4 Bergman kernel 

We use the notations of [H], and we refer the reader to [H] for a more detailed description. 

On every compact Riemann surface £, with a given symplectic basis of non trivial cycles , is defined uniquely 
a 2nd kind differential called the Bergman kernel [2] B(pi,p2) (which we regard as a 2nd kind differential in the 
variable p\ S C), that satisfies 

i) B(p\,p2) has a double pole, with no residue, when p\ — > and normalized such that 

dx(pi)dx(p 2 
'pi— P2 (x(pi) - x(p 2 )) 

where x(j>) can be any local parameter in the vicinity of p 2 . 

ii) 

B( Pl ,p 2 ) = 



B(PX,P2)_ ~_. ,17,^ +^ ( 4 - 13 ) 
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It is easy to see that the Bergman kernel is unique, because the difference of 2 Bergman kernels would have 
no pole and vanishing A— cycle integrals, i.e. it would vanish. 
More explicitly we have: 

B(px,P2) = d Pl d P2 log(0(u(pi) - u(p 2 ) - «)) (4.14) 

where 8 is the theta-function, u(p) is the Abel map, and k is some odd characteristics. 

For example, in the case the spectral curve has genus zero (the so-called 1-cut case), C is the Riemann sphere, 
i.e. the complex plane with a point at oo, and B(pi,p 2 ) is the meromorphic bilinear form B(pi,p 2 ) = fof l^sa ■ 
Another example is the case where £ is a torus of modulus t: C = C/(Z + rZ), for which the Bergman kernel 
is the Weierstrass function: B(pi,p 2 ) = (p(pi —p 2 ;r) + C)dp\dp 2 . 

4.5 Third kind differential 

For any p € C and two points q and o in the fundamental domain, we define: 

dS q , (p) = [ B(p,q') (4.15) 

J o 

where the integration contour lies in the fundamental domain (i.e. it does not intersect any A or B cycles). 
dS q>0 {p) is a meromorphic differential form in the variable p, whereas it is a scalar function of q and o. It has a 
simple pole at p = q with residue +1 and a simple pole at p — o with residue —1: 

Res dS q (p) = +1 , Res dS q (p) = -1 (4.16) 

p^q ' p^o 

i.e. it behaves locally like ;^gj~7^ when p — > q, in any local parameter x(p). Moreover it has vanishing A cycle 
integrals: 

j dS q , o = (4.17) 

Since it has only one simple pole in the variable q, this 3rd kind differential is very useful for writing Cauchy 
residue formula. For any meromorphic differential form u>(p) we have: 

uj{p) = - Res dS q , (p) Lu(q) (4.18) 

q~tp 

and, using Riemann bilinear identity |25U2G|. if V«, §. u> = 0, and u> has poles a^'s, we may move the integration 
contour and get: 

u(p) = V Res dS q . (p)Lu{q) (4.19) 

* — * q—>ai 

i 

This identity was the main ingredient in solving loop equations for the 1-matrix model in |16| . 



5 Solution of the Loop Equation 

In this section we solve the loop equation to all orders in the topological T 2 /N 2 expansion. 

We first need a technical lemma which consists in proving that the solution is unique, and then we use this 
uniqueness to try a guess similar to that introduced in [9] which makes the loop equations easier to solve. 

We find the one point resolvent and the k point resolvent for the first matrix of the chain, and in fact we 
find that they coincide with the correlators defined in [19] for the spectral curve 
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5.1 Unicity of the solution 

Equation (|3.10fl fixes the large N/T expansion of Wo(^i(p)m Take equation (|3.10j) and substitute the 



expansion of W 2 .i(xi,x[), Wq{x-\), U(xi,x 2 ) and E(xi,x 2 ). Then to order jt^k we obtain 



JV 2 

h 



(c h2 x 2 ~Y(x 1 ))m h \x 1 ,x 2 )+W^\x 1 )m o \x 1 ,x 2 ) = 

= E (h \x u x 2 ) - W^- m \x 1 )U^\x l ,x 2 ) - W^- 1) (x 1 ,x 2 ;x 1 ). 

Suppose you know U^ h ^(x± 3 X2 ), W^ h '\xi) and E^'\ Xl ,x 2 ) for h! < h. We prove that we can find those three 
functions for h! = h. Consider x\ = x±(q) and x 2 = x 2 (p) (and so Y{x\{q)) = ci t2 x 2 (q)) with p and q living on 
the algebraic curve 

ci, 2 (x-2 (p) - x 2 (q)) m ( Xl (q),x 2 (p)) + W^ h) (x, (q))U^ (a* (g), x 2 (p)) = 

= m{x 1 {q),x 2 {p)) - £ W^ h - m \x 1 (q))m m \x 1 (q),x 2 (p))-W^- 1 \x 1 (q),x 2 (py,x 1 (q)). ^ 

m— 1 

Begin with h = 0. Consider the solutions for the equations E(xi(q),x 2 (p)) = 0. For every there are 

D 2 + 1 different solutions y(xi(gW)) (sitting at points that we call g^,^ 1 ),...,^) on the curve, with the 
convention that q^ = q). Then we can write 



D-, 



m ( Xl (q), x 2 (p)) = K ft (c h2 x 2 (p) ~ Y( Xl (g«)) 

i=0 



(5.3) 



i=l 



where the constant if is derived in the next section. Recall that Xi(q^) = xi(q^) but in generaQ Y(xi(q^)) ^ 
Y(xi(q^)) for i ^ j. 

Consider now h > 0. 
Write now the equation for an arbitrary h and take p — > g° = g 

^ W (x 1 (q))L>( )(^i(<z),x 2 (g))= J B(' i )(x 1 ( (Z ) I x 2 (q)) 

- E ^"^Ms))^" (11(g), a; 2 (9)) - wfr'W^^g);^))- 

m— 1 

This equation shows (by recursion) that W"d^(xi(g)) is a meromorphic function on the spectral curve, and 
because of our hypothesis, it has poles only at branch-points, and it has vanishing A cycle integrals. Let us 
write Cauchy residue formula (|4.18j) : 

W^ix^dx^q) = - Res dS q> . {q) W^ix^dx^q') (5.5) 

q'~>q 

Using Riemann bilinear identity, and the fact that both dS and W^dx have vanishing A cycle integrals, we 
can move the integration contour and get <|4. 19[) : 

W^ h \ Xl (q))d Xl (q) = V Res dS g ,, (q) W^ h) {x 1 {q'))dx l {q') (5.6) 



3 As we show later all the fc-functions of the type 
W. 



i*u (»iw. ■ ■ • . = (n 9V/ J (g(i))) ) wb (h) (*i w) 

can be determined from the equation l|3.10ll exactly in the same way as in |19ll9]. . . . 

7 The function Y(xi(q^)) is multi valued in the x\ plane. On the other side on the algebraic curve it is not multi valued. The 
index i indicates precisely different xi-sheets, and thus different values of Y(x%). 
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Now we replace Wq (xi(q')) in the RHS with the loop equation l5,4} and using that E w (x 1 ,x 2 ) is a polynomial 
and has no poles at finite x\, that U^\xi{q), x 2 {q)) vanishes at most as a square root at the branch points and 
that Ax\(p) vanishes linearly at the branchpoints, we find: 



W^ h \x 1 {q))dx 1 {q) = V Res \x 1 {q'))dx 1 {q')AS ql >0 (q) 

* — ' o'— >a 



q — >oc 

h-1 



Z ^^^UW{x 1 {q'),x 2 {q'))\^ 1 ° ^ 7 > 

+ W^- 1 \x 1 (q'),x 2 (q'y,x 1 (q')yj. 

where everything on the RHS is known from the recursion hypothesis, and thus determine uniquely WQ h \xi(q)). 

Then, consider again equation (|5.4p and find E^ (x 1 (q), x 2 (q)) (equal to E^ (xi(q^) , x 2 (q)) by the defini- 
tion of gW) 



m{ Xl (qV),X 2 {q)) = E^(x 1 (q),x 2 (q)) = {x^U^ {x x {q),X 2 {q)) 



h-1 (5 g\ 

- E W^ h - m \x 1 (q))U( m \x 1 (q),x 2 (q))-W^- 1 \x 1 (q) > x 2 (qy ) x 1 (q)) 

m—1 

and reconstruct E^(xi(q),X2(p)) using the Lagrange interpolation formula 

E w (xi(q),x 2 (p)) = y — — - J -£ , . 5.9) 

i Uj^iiMQ ) -x 2 (q w )) 

Finally equation l|5.2p gives U^ h '(xx(q), x 2 (p)). 

Therefore we have proved our recursion hypothesis to order h. 

All this procedure allows us to solve recursively the master loop equation, thus indicating that the solution 
is unique once EW(x u x 2 ) (or equivalently Y(xi(p)) and x\{p)) is given. We could iterate this procedure 
indefinitely. We now show a much better way to solve the master loop equation. 

5.2 Solution of the equation 

The solution being unique, we only have to find one solution. The equation 

rp 2 

jj 2 -W 2; i(xi,x 2 ;xi)+(cx t2 x 2 - V{{xi) + W {x 1 ))U(x 1 ,x 2 ) = 

= - Wi{x\,x 2 ) + {V((x±) - ci, 2 x 2 )S(xi,x 2 ) = E(xx,x 2 ) 

is indeed solved by the expressions 

/D 2 f 

E( Xl (p),x 2 ) = -K» ( [] c i,2-^2 - V({x x {p)) + -tr 



u=o x 

/D 2 



N \x±(pW)-M 



(5.11) 



U( Xl (p),x 2 ) = -K» ( [] ( c Xj2 x 2 - V;( Xl (p)) + ^tr 



N \xi(pW)-M 



and can be proved following [9]. The product runs over the D 2 + 1 sheets of the algebraic curve viewed from 
the x\ variable point of view. The Oth sheet is by definition the sheet in which the point p is sitting (that is, 
p = p( '). The notation " (•••)" means that if we expand the product in cumulants, the connected two point 
correlators must be replaced by W -i(xi(p); xi(p')) = W -i{xi(p);xi(p')) + ( Xl ( p )- Xl ( p >)^ ■ 

These expressions are not of practical immediate use, but if we expand them in powers of x 2 they reveal the 
equation that lead us to the explicit solution. All the information is contained in the highest powers. 
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(C1,2X 2 ) 



D 2 + l 



1 

Cl,2 



d 2 ---d„5 



(2) \ d 3 ---d n $ . /g) v d4— d„s 

f 5<j 3 +i \ 



C 2 ,3 



3d 2 + l 



V C3 - 4 J 



Cn,n+1 



= if 



(5.12) 



(ci,2Z2)' D2 : 



A' 



Vi(xi) -d 3 ---d n s 



(2) 

C 1.2.9rf 2 

o (2) 



D-2 



i=0 



N \xi(pW)-M 



(5.13) 



(2) 

■9rf 2 + l i=0 



= X)y(si(p W )) 



where we have defined as usual V/(xi) - W ( Xl (p^)) = Y(xi(p®)) 

{Cl,2X2) ■ 



P(xi) -d 3 ---d n s 



(2) 

c l,2grf a 

(2) 
9d 2 + l 



v;( Xl ) - (d 3 ■ ■ ■ d n5 - 1) 



(2) (2) ' 

1 c h2 g d2 ' Ci, 2 fli 



2 (2) 

yd 2 +i 



(2) 

5rf 2 



r>2 



igfr(x 1 (/ i )))F( a;i (^)) + ^^o;i(x 1 (^ 



(5.14) 



where P(xi) = PoL^ y/(a;i)Wo(.Ti) was already defined in <|2. 14|) and 

Wo ; i(n(p),a;i(g)) = Wb;i(xi(p),xi(g)) + 



1 



(xi(p) -xi(q)) 2 ' 

refers to the substitution mentioned above for two point correlators. 



(5.15) 



The equation (|5- 12[) allows us to determine the constant K. The equation l|5.13p allows us to modify the last 
equation. When doing the T 2 /N 2 expansion, equation (|5 . 13[) implies 



D 2 



W W (zi(p w )) = for h > 



(5.16) 



i=0 



Apply 



dVi{x{q)) 



to equation (|5 . 13jl we find also 



D, 



1 



i=0 



(xi(g) -xi(p)Y 



(5.17) 



Using all these equations we find that equation (|5.14p can be transformed into 



D 2 

E 

i=0 



Y( Xl (p«)) a + ^Wo-Ax! (p® ) , Xl (pW )) 



= (^'(.tx)) 2 - P(x x ) - <fe . . . d n S (ds . . . d„5 - 2) 



' C l,29d 2 
(2) 

> 3d 2 +i , 



2^3 • • • d„s 



. 2_(2) 
c l,2 gja-l 

(2) 
9d 2 + l 



(5.18) 



Expanding the equation in j£ = ft as in [9] we get for /i > 1 the equation (with y(p) = Y(xi(p))) 



D 2 



2Y,y(p {l) )w^ h \x 1 (p^)) 



1=1 



D 2 h-1 



D 2 



= E E w-o (m) (»i(p (i) ))Wo (h ~ ro) («i(p c<) )) +E<r 1) (^(f w )' a; i(f w )) +2P (fe) (xi) 



(5.19) 
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The rest follows exactly the same lines as in [9]. We will however recall the main steps. Let us define the 
following meromorphic differentials from the correlation functions 



d 



4»(p u ..., Pk )= dx 1 ip j )j (n dVlixliPl)) ) ^Wi)) (5.20) 

and rewrite equation (|5.19[) as 

D 2 

2^y(pW)^ ) (pW)dx 1 (p«) = 

i=l 

= E E ^r\p {i) )^ m \p {i) ) +E4 M (p (i) .p (i) ) + 2P^)(x 1 ( P ))d a;i b) 2 



(5.21) 



Define also the third kind differential dE Pt p(q) = dS p , (q) — dSp t o(q), where p is the conjugated point of p. 
Finally apply the operator ^ Q Res p ^ Q I y ( Xl (pjf-yfx! (p)) ( wnere a are the branch points of the curve) to the 
equation l|5.20p . After some algebra we find 



'h-l 



= -£ Res I rj^f^ M E 4 m \p)"i h - m) (P) +4 h - 1) (p,P)) (5.22) 
^P^a 2 (y(p) - y(p))da;i(p) > 



V?n — 1 



which is the first of a tower of recursion relations. The rest can be obtained by applying the loop insertion 
operator to this first one and reads 



(5.23) 

m=0 JC-R" / 



where {pk} is a collective notation for k points on the curve, and K = {1, . . . , k} is the set of indices. In the 
expression, J stands for a subset of j elements of K, K\J for the complement of J in if and the sum over J 
and to counts all different subsets and genus, except ( J, m) = (0,0) and (J, to) = (K,h). 

Therefore we have found that the meromorphic differentials (Si) • • • > Sfe) satisfy exactly the same recursion 
structure as those of [19], thus all manipulations done in [19] and other references therein that depend only on 
this resursion structure need not be repeated here and can be taken as a fact. For instance it was shown in [19] 
that the recursion for the F g 's and W^'s can be represented digrammaticaly, and so the same happens here. 

6 Moduli of the chain of matrices and topological expansion of the 
free energy 

In order to find the free energy it is important to understand which are the moduli of the chain of matrices, 
and how they change when we change the curve (always within the matrix chain moduli space). 



6.1 Moduli of the chain of matrices 

The chain of matrices is completely characterized by the potentials Vi(x), . . . , V n (x), the interaction parameters 
Ci,i+J5, the temperature parameter T, the eigenvalues and multiplicities of A and the filling fractions a. 

8 Note that c n , n +i can be absorbed in a redefinition of A. We will see later how this appears in the moduli variations. 
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It is clear that we can express all these parameters in terms of the meromorphic functions on the curve 
Xi(p), . . . ,x n (p) as follows (with notations borrowed from [4} [191 15]): 

£i =tt~ <£ Ci, 2 x 2 dxi = (f c k . k+1 x k+1 dx k = -— !- / c k . k+1 x k dx k+1 
2tt« J a, 2m JAi 2m JAi 

too =T = Resci ;2 x 2 dxi = Res c k , k+ ix k+ i&x k = -Resc k , k+ ix k dx k+ i 

oo oo oo 

h 

U = - —T = -Resc„,„ + ix n dx„ + i = - Rcsc k . k+ ix k dx k+1 

= Res c ktk+1 x k+ idx k = Res c lj2 x 2 dxi (6.1) 
gy =ciq Resx^x 2 da;i 

J oo 

j > 2 : g) k > =c k , k+1 Resx^.Tfc+idxfc = c fe _i k Res x7 J a; fc _idx fe Vi 

00 A s 

tt.Ai = -t?.x n+ i(Aj) = -c„. n+ iResx n+ i(p)x n (p)dx„+i(p) 

Ai 

Something deserves attention here: note that can be expressed in s + 1 different ways by changing which 
pole Ai or 00 we consider. As we will see later, in order to stay within the matrix chain moduli space, any 
variation of the curve around one of these points should bring associated other variations around the other 
points so that the new g's can still be obtained from any of them. Also note that we have not specified how 
to obtain Cj^+i. They appear in the other equations as to indicate that they are free to choose. Indeed these 
parameters can always be absorbed into the other parameters of the model (as the equations above indicate) . 
It can also be viewed as a rescaling of the meromorphic functions x% (p) . 

Now, we study how the spectral curve changes when we change these parameters (or vice versa). 
Let us define the variations fl of the curve by their effect on the differential 01,2X2 (p)dxi(p). Variations of func- 
tions or forms, are defined with respect to some fixed variable. There is a Poisson-like structure (thermodynamic 
identity) indicating how to relate variations with respect to different fixed parameters. The meromorphic form 
O is defined as: 

<Mci,2X 2 (p)dxi(p))| xi(p) = Mci, 2 x 2 (p))| p dxi(p) - fo(xi(p))| p ci )2 dx2(p) = -Q(p) (6.2) 
In general we want fl to be written in the form 



n(p) = / B(p,q)A(q) (6.3) 
where dfl is a path which does not intersect circles around branch points. 

6.2 Variation of filling fractions 

For variations of the filling fractions we choose 

fi(p) = -2iTrduj(p) = - I B(p,q) (6.4) 

JBj 

so that d£l = Bj and A(q) = — 1. From (|6.ip we have 

Snei = 5 jt , 5 n ta = , 5 n g^ n) = , 8 n \ t = (6.5) 
so that indeed, 5-2wrdiy = ^f~- Using Theorem 5.1 in [l9j we can write 

-^w k h) (pi,...,p k ) = -j( w ( k %(pi,...,p k ,q) (6.6) 
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6.3 Variation of the temperatures 

Similarly we define for = T and . , 



flip) = -dS a , a , = / B(p, q), i.e. SO = [a, a'],A=l 

J a 



3.7) 



where a, a' € {oo, Ai, . . . , A 5 }. This variation produces the following modifications of parameters 

8 n ei = , Satff = 5 a!0 - 5 a > t p , 5ng { - n) = , 5 n X t = (6.8) 

which can be written as <5-dS a a , = gf g§-. This makes sense since ^ Q t a = 0. Again Theorem 5.1 in [19] 

enables us to write 



d 



W a ~ W) W k' ) ^^--'Pk) = j «>fc+i(Pl»" (6-9) 

6.4 Variation of the potentials 

Observe that if we don't consider variations in c^i+i we have 

ci, 2 (5x2-dxi - 5x\.dx 2 ) = 5x 2 .(V 2 (x 2 )dx 2 - c-z^dxz) - {(8V 2 ){x 2 ) -dx 2 + V 2 (x 2 )6x 2 .dx 2 - c 2 , 3 8x 3 .dx 2 ) 
= -(5V 2 )(x 2 ).dx 2 + c 2 ^{8x 3 .dx 2 — 5x 2 .dx 3 ) 

(6.10) 

= -{SV 2 ){x 2 ).dx 2 (SV^)(x n ).dx n + c n ,n+i(Sx n+1 .dx n - Sx n .dx n+1 ) 

In particular, if the Ai are kept fixed, the last term 5x n +i-dx n — Sx n -dx n +i has no pole at Aj. 
Variations of V\ 

If we vary only V\, more precisely if we vary only gf\ we see that O = ci t2 (6xi.dx 2 — 5x 2 .dx\) has no pole at 
the A,'s, and near oo, if we work at fixed x±, we have c\. 2 5x 2 ~ $9j x{ 1 + 0(x^ 2 ), and in addition, since we 
don't vary the filling fractions, we know that §, O = 0. All these considerations imply that the variations of 
V\(x) are given by the same formulas as in [9]: 

0(p) = -Boo,j(p) = - Res B{ Pl q)x\{q) (6.11) 

J q — >oo 

thus 50 is a small circle around 00 and A(q) = ^ Xl ^ ■ With this variation it is easy to check that 

8nei = , S n t a = , 6 a g£> = d k ,iSj, m , S n X t = (6.12) 
and so we can say that (5_b , = fn , and from Theorem 5.1 in 1191 : 

9 ('')/ \ 13 Kite) 3 (h) , S ,r. 10 N 

—7yrw k '(pi,...,Pk) =Res ; — wW 1 (pi,...,Pk,q). (6.13) 



dg 

Variations of V 2 , . . . , V n 

For the other potentials Vk with 2 < k < n, if we vary g\ k \ near 00, at fixed xi, we have Sx 2 = ci t2 S(V{(xi) — 
T/xi + 0(x]~ 2 )) = 0(x^ 2 ), therefore O has no pole at 00. We have seen that the pole of O at Ai is given by 
6(y£)(xk)dxk, therefore near Ai we have O ~ —x^T^dxk. This implies that 

O(p) = -J2 B k k j(P) = -~ E Res B ^ 9)4(9) ( 6 - 14 ) 
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thus dn is a contour which surrounds all A.; (and no other poles), and A(q) = —-^ ^j ■ Then, we can say 
that <5_ v b- — — iTT i and from Theorem 5.1 in 1191 : 

— 1ja w I H \pi> ■■■,Pl) = VRes w^ipi, ...,pi,q). (6.15) 

6.5 Variation of the Aj's 

Similarly, we see that when we vary Ai, fl has no pole at oo, and near Aj, it behaves like —dx n . Therefore we 
have: 

fi(p) B X » = Res Bfog^fa) =2^ (6.16) 

JV q^A; iV dx„+l(Ai) 

-^-w i j h) (j)i,...,pi)^T^ Rpsw$ 1 (p 1 ,...,pi,q)x n (q). (6.17) 

6.6 Variation of the c^k+i 

Again, allowing variations in Ck.k+i we find that ft behaves like Xkd.Xk+i, therefore 

fi(p) = = l]ResS(p,g)^( g )^ +1 (g) (6.18) 



d 



~i (pi,---,Pi) = *52 R esw i+i(Pi> ■ ■ ■ ,PhQ) x k(Q)xk+i(q) • (6.19) 

OCk,k+l ~ Xi 

6.7 Summary of moduli 

Using Cauchy formula, we may write: 

ci,2X2(p)dxi(p) = -Res dS q , (p) ci ox 2 (q)dxi(q) (6.20) 

Then, we move the integration contour, and we take into account the boundary terms using Riemann bilinear 
identity, we get: 

ci,2X2ip)dxi(p) = Res dS q ^ {p)ci^X2{q)dxi{q) + 2i-K'^eidu i (p) (6.21) 

q — >oc,Ai 

l 

The residues near the poles oo and near the Xi are computed by the local behaviors. 
• near oo, we have X2 ~ V{{x\) — jj- + 0(x^ 2 ), i.e. 

Res dS q , {j)) c 1 oX2{q)dx\(q) = Res dS„ (p) dVi(x\{q)) - T Res dS q , {p) dxi ^ 

q^oa q—*oo q—>oo Xl{q) 

= - Res B(q,p)V 1 (x 1 (q)) +TdS 00 , (p) 

q — >oo 

(1) 

= -E — Res B ^p) (xi( q )y +Tds 0C . (p) (6,22) 

3 

= E -— B °o,j(p) + TdSoo^ip) 
j 3 
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• near Xi, we have 

ci a x 2 dxi = Ci t 2(d(xix 2 ) - x 1 dx 2 ) = ci t2 d(x 1 x 2 ) - V 2 '(x 2 )dx 2 + c 2t3 x 3 dx 2 



(6.23) 



= d(a i2 xix 2 - V 2 (x 2 ) + c 2t3 x 2 x 3 - V 3 (x 3 ) . . . 

~t~ Cn—l,nXn—lXn Vn(x n *j -\- C n ^ n ^-iX n X n ^-i ) C n ^ n +iX n dx n -\-i 

~ d{ci }2 xix 2 - V 2 (x 2 ) + c 2 , 3 x 2 x 3 - V 3 (x 3 ) . . . 

+ c n -i. n x n -ix n - V n (x n ) + c n , n+1 x n x n+ i) - T-^- dXn+1 + £>(() 1) 

N x n+ i - Xi 

Therefore we have: 

ResdSq. {p) c lt2 x 2 (q)dx 1 (q) = 

= ~ T Tf dS \ o(p) ~ Res dS qAp) d{V 2 H \-V n - ci t2 xix 2 c n . n+ ix n x n+ i) 

= dS \ o(P) + Rcs B (^P) ( V 2 H \-V n - ci. 2X1X2 c„.„ + ia;„x„ + i) 

, n (k) n 

= ~T^dS x . » • Rc - S B ^ ( Xk bW ~ E C W Re - S S( g ,p)x fe ( g ) a;fc+1 ( g ) 

fc=2 3 J <^ A ' fe=l «^ X< 

All this can be summarized as: 

k=2 j 1 



Tl n_1 

i k.—l i 



(6.25) 



Notice that c„ jra +i does not appear, in fact the term that would logically give the associated contribution, 
it is better used to encode the variations of Xi. It is clear that the Xi contain already the information of Cn, n +i- 

6.8 Topological expansion of the free energy 

With all that information we are now ready to derive the free energy The free energy In Z = T = 

J2 g (N /T) 2 ~ 2g 1F^ 9 \ is determined by its derivatives: 

- <Tr(M k y>=--^. (6.26) 

The result that we wish to prove is that: 

T^=F g (E^) (6.27) 
where F ff 's are the symplectic invariants of [19], for the spectral curve E^ ' . In particular for g > 2 we have: 

Fg(E {0) ) = ~z — ~tt" Y Resw^ $ , d$ = c h2 x 2 d Xl (6.28) 
z — Zq ' — ' a 

The expressions of Fq and Fi are a little bit more difficult to write [3H [31], and we refer the reader to [19j. 
Notice that when there is no external field, i.e. A = 0, was already computed in [15], and it coincides with 
F . 

The Fg's of [19] have the property, that under any variation ft, we have: 

5F g = f w[ 3 \q)A(q) (6.29) 
Jon 
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In particular with k = 1, it proves that 

dF g 

,(D j ^ g g (D 



dF a i M , aj^ff) 

- Res w\ 9) x{ = — -— (6.30) 



Then, we prove it by recursion on the length of the chain n. 

The n = 1 case was done in [19]. Now, assume that it is true for n — 1. 

We have just seen that J-^ — F g is independent of V\, therefore we may compute it for the case where V\ 
is quadratic. When V± is quadratic, the integral over the first matrix of the chain, Mi, is a gaussian integral, 
and Mi can be integrated out, so that when V± is quadratic we are left with a chain of n — 1 matrices, and 
we get Tn = Fn-i- From the recursion hypothesis, we have T^l\ = F g (E^ {x2, £3)), and one should notice 
that the F g 's of [19] have the symplectic invariance property, i.e. they are unchanged if we make a symplectic 
transformation of the spectral curve, or in other words, if we add an exact differential to ci, 2X2(&Ei . In particular 
we may work with 02,3X3^X2, and thus F g (E^ (xi, X2)) = F g (E^ (x2, X3)). This proves the result. 



7 Other Considerations 

In the previous two sections, we have solved the loop equations to all orders, and we have found that the solution 
is given by the symplectic invariants introduced in [19], for the spectral curve (xi,x 2 ). As a consequence, 
all the properties studied in [19] apply. 



7.1 Symplectic transformations 

Remember that the spectral curve E(xi,X2) = is equivalently given by the data of two meromorphic func- 
tions xi(p), X2{p) on C. Indeed, given two meromorphic functions, it is always possible to find a polynomial 
relationship between them. We shall write the spectral curve: 

E h2 = {(x 1 (p),X2(p)) / peC] = (xi,x 2 ) (7.1) 

Since any Xi is a meromorphic function, we can also define the following algebraic spectral curves: 

Ei,j = {{xi(p),Xj(p)) j peC} = (x^Xj) (7.2) 



It was found in |19[I20|. that the F g 's are unchanged under symplectic transformations of the spectral curve, 
for instance if we add to xi any rational function of X2, or if we exchange xi <-> X2, or if we change xi — > —x\. 

For instance we could change c\.2.x\ — > 01,2X1 — V{(x2) = —02,3X3, and then X3 — > —X3, and then recursively 
d,i+iXi — > Ci,i+ix ? ; - V{ +l {xi) = -Ci+i : i+2Xi+2- This shows that: 

F g = F g (E i>i+1 ) = F g {E i+lti ) Vl<i<n (7.3) 

However, one should keep in mind that the correlation functions are not conserved under symplectic trans- 
formations, only the F s 's are. 



7.2 Double scaling limits 

We have seen that as long as the spectral curve is regular (all branch-points are simple), the F g 's and all 
correlation functions can be computed, and it was found in [l9j that they diverge when the curve becomes 
singular. This type singularities were already found in the one and two matrix model, and in [19] for generic 
spectral curves, but it is still important to show that it appears in this context too. 

It was found in [19], that if the spectral curve depends on some coupling constant (T for instance), if the 
spectral curve develops a cusp singularity at say T — T c of the form 

y ~ x p/q (7.4) 
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then the F g 7 s diverge as 




of the singularity, and which is the spectral curve of the (p, q) minimal model jl ll [29] . All this is detailed in 
|19j and we refer the reader to that article for more details. 

As usual, singularities of formal series are related to the large order asymptotic expansion of the general term 
of the series [33] , and the double scaling limit is thus related to the asymptotic enumeration of large discrete 
surfaces, and in some sense to their continuous limit, i.e. Riemann surfaces. A (p, q) minimal model may occur 
as soon as two of the Vi's have degree larger than p and q [11]. Here, we see that the double scaling of the chain 
of matrices describes a (p, q) minimal model on a random lattice. This is related to the Liouville conformal field 
theory coupled to minimal models (p, q) . This phenomenon is expected [30], HH [321 E] an d is already known to 
be present in the one and two matrix models and more generaly in [19j. 

7.3 Modular transformations and holomorphic anomaly equations 

In order to compute the F g 's and the solution of loop equations, we have made a choice of cycles A, related to 
the choice of the minimum around which the formal matrix integral is defined. However, it is interesting to see 
what happens if one makes a different choice of cycles, i.e. if one makes a modular transformation. This was 
studied in [19] and [22]. 

A modular transformation changes the Bergman kernel B(p, q) with a constant symmetric matrix k: 



(where n t j — § B dui is the Riemann matrix of periods of C) then the Bergman kernel is called Schiffer kernel 
and is modular invariant. 

More generally, the modular transformations were computed in [19] . and they satisfy the so-called holomor- 
phic anomaly equations, and that gives a strong support to the Dijkgraaf-Vafa conjecture that matrix models 
are topological type B string theory partition functions [3] [22]. 

7.4 Convergent matrix integrals and filling fractions 

So far, we have considered formal matrix integrals, defined by expanding the integrand in the matrix integral, 
near a given extrema specified by a set of filling fractions. We worked at fixed filling fractions. 

On the other hand, convergent matrix integrals should correspond to integrals over (Hn)" 1 . The integration 
path can always be written as a linear combination of steepest descent paths (those used for formal integrals) , 
and the full convergent matrix integral is obtained as a linear combination of formal matrix integrals. More 
precisely, the convergent matrix integral should be a sum over filling fractions of the formal ones. 

The summation over filling fractions was computed in [21], and just amounts to multiplication of the formal 
matrix integrals by a theta function. We refer the reader to [6] [21] for more details. 

8 Limit of a continuous chain of matrices 

In this section, we briefly explore some consequences of our method for the continuous chain of matrices. 

The "matrix-model quantum mechanics", is obtained [14] as the limit n — > oo, and with the choice: c^j+i = 




(7.6) 



where dui are the holomorphic forms on C such that §. dui = Sij. 
In particular, if 




(7.7) 



-, and: 



1 

Ci,?:+i — - 



e 




Vi(x)=eV(x,ei) + ^ 



(8.1) 
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The index i is rescaled as a continuous time t = ei: 

t = ei , 0<t<t f = en (8.2) 

The matrix integral thus becomes: 

Z= [ D[M(t)] e-*/o / * Tr[V(M(t),t) + iM(tf] (8 _ 3) 



The spectral curve is determined by the equations V-(xi) = Ci,i+iXi+i + dj-xXi-i which become Newton's 
equation of motion [14] to leading order in e: 

V'(x,t)=x{t) (8.4) 

and the resolvent of the first matrix is: 

W{x, 0) = V({X!) - c h2 x 2 ~ -x(0) (8.5) 

The topological expansion is thus: 

Z = e^™"' F » (8.6) 

F g = F g (£(t)) , £(t) = (x(t) 7 -x(t)) (8.7) 

The spectral curve £(t) = {x(t), —x(t)) is thus the dispersion relation, i.e. the relationship between velocity 
and position, it may depend on the time t, but from symplectic invariance, we see that F g (£(t)) is a conserved 
quantity, independent of the time t. 

For example, if the potential V(x,t) = V(a;) is independent of t, the kinetic energy K is conserved and the 
dispersion relation is: 



\x 2 - V{x{t)) = K 



and the spectral curve is: 



£(t) = (x(t), y/2(V(x(t)) + K)) (8.9) 



Consequences of those relations need to be further explored, and we leave the continuous chain of matrices 
for another work. 

9 Conclusion 

We have computed explicitly the topological expansion of the chain of matrices with an external field, and we 
have found that the F g 's are precisely the symplectic invariants of |19j . 

We have also computed some of the correlation functions, but not all of them, in particular we have not 
computed mixed traces (which count discrete surfaces with non-trivial boundary conditions) . Mixed traces were 
computed in the 2-matrix model case in [20l [17] , and it would be interesting to see how that could be extended 
to the chain of matrices. 

We have also briefly started to explore the limit of matrix quantum mechanics, i.e. the limit of an infinite 
chain of matrices, but this topic needs to be studied in deeper details. 
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